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Q ■ Abstract 

The relationship between the properties of nuclear matter and structure 
functions measured in lepton-nucleus deep inelastic scattering is investigated 
using light front dynamics. We find that relativistic mean field models such as 
the Walecka, Zimanyi-Moszkowski (and point-coupling versions of the same) 
and Rusnak-Furnstahl models contain essentially no binding effect, in accord 
P^ with an earlier calculation by Birse. These models are found to obey the 

Hugenholtz-van Hove theorem, which is applicable if nucleons are the only 
degrees of freedom. Any model in which the entire Fock space wave function 
can be represented in terms of free nucleons must obey this theorem, which 
implies that all of the plus momentum is carried by nucleons, and therefore 
that there will be essentially no binding effect. The explicit presence of nu- 
clear mesons allows one to obtain a modified form of the Hugenholtz-van Hove 
theorem, which is equivalent to the often-used momentum sum rule. These 
results argue in favor of a conclusion that the depletion of the deep inelas- 
tic structure function observed in the valence quark regime is due to some 
interesting effect involving dynamics beyond the conventional nucleon-meson 
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treatment of nuclear physics, 
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I. INTRODUCTION 

The European Muon Collaboration (EMC) effect in which the structure function of a 
nucleus, measured in deep inelastic scattering at values of Bjorken x > 0.4 corresponding 
to the valence quark regime, was found to be reduced compared with that of a free nucleon 
was discovered almost twenty years ago |1|]. Despite much experimental and theoretical 
progress ||||, no unique and universally accepted explanation of the depletion has emerged. 
The immediate parton model interpretation that the nucleon bound in a nucleus carries less 
momentum than in free space seems uncontested, but determining the underlying origin 
remains an elusive goal. 

One popular explanation is that conventional nuclear binding effects are responsible. The 
conventional lore is that the nuclear structure function F 2 a{x) (which gives the momentum 



distribution of a quark in a nucleus as a function of the fractional momentum carried) can 
be obtained from the light front distribution function f(y) (which gives the probability that 
a nucleon carries a fractional momentum y) and the nucleon structure function of a free 
nucleon F 2N using the relation M: 
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dyf(y)F 2N (x/y). (1.1) 



This formula has a simple interpretation as an expression which gives a probability as a 
sum of products of probabilities. The variable x is the Bjorken variable x = Q 2 /2Mi>, 
and y is the A times the fraction of the nuclear plus-momentum carried by the nucleon, 
The plus component of a four vector is the sum of the time and third spatial component, 
so if k^ is the momentum of a nucleon and P M is the momentum of the target nucleus 
y = (k° + k 3 )A/P + = (k° + k 3 )A/M A = (k° + k 3 )/M, in which the nucleus is taken to 
be at rest with P + = Ma- One can easily use conventional nuclear physics to obtain 
the probability that a nucleon carries a three momentum k, but, if one uses only naive 
considerations, one faces a puzzle when deciding how to choose the value of k°. Should one 
use the average separation energy, or the average nucleon mass M, or possibly the effective 
mass in the chosen many-body theory? 

The essence of the binding explanation is that k° is given by the free nucleon mass M 
minus the average separation energy e. Then f(y) is narrowly peaked at y — 1 — e/M (e 
(with ~ 70 MeV for infinite nuclear matter ||). In this case, the structure function of a 
bound nucleon is approximately obtained by replacing F 2 n(x) by F 2 n{x/(1 — e/M)). The 
increase in the argument leads to a significant reduction in the value of the nuclear structure 
function The theoretical understanding of the binding effect (as of 1996) is reviewed nicely 
in the book |3|, which summarizes the various treatments as "not completely satisfactory". 
This kind of explanation seems very natural because nuclear binding is known to occur, so 
such an effect must be understood thoroughly before hoping to extract information about a 
possible host of more interesting exotic effects. In any case, one needs to supply a derivation 
to avoid the need to arbitrarily choose a prescription for k°. 

This need drove one of us on to the light front ]7|,|8| . That is, to attempt to use light front 
dynamics to derive the nuclear wave function. The reason for this is that, in the parton 
model x is the ratio of the plus component of the momentum of the struck quark to that 
of the target, and it is the plus component of the momentum which was observed to be 
depleted by the EMC. In the view of Ref. ||, using light front dynamics is the most effective 
way to assess the influence of binding effects. However, one must pay the price of computing 
nuclear wave functions using these dynamics. 

The first attempts 0Q in this direction employed the popular and successful Walecka 
model 0] which has many effective descendants fTTHTBl. The salient result was that vector 
mesons carried 35% of the nuclear plus momentum and nucleons only 65% (P^/P + = 0.65), 
far smaller than the value (1 — e/M) ~ 0.95 needed to reproduce the observations for the iron 
nucleus. However, the connection between the nucleon momentum distribution computed 
using light front dynamics and that used in computing the deep inelastic structure function 
was not made. Recently, the authors of Ref. WM have claimed that quark distribution 



functions are not parton probabilities. Their message to us is that, in any situation, one needs 
to derive the connection between the constituent distribution function and the observed 



data. That work stimulated us to undertake the present investigation in which we derive 
the connection between the nucleon momentum distribution and the structure function 
measured in deep inelastic scattering. 

First we outline our procedure. We start in Section [TJ by considering relativistic models of 
infinite nuclear matter computed using the mean field approximation. We derive and apply 
the nucleon distribution function /jv(|/) appropriate for use in computing deep inelastic 
scattering structure functions. The function /jv(y) is shown to be the one which maintains 
the covariance of the formalism, and in which the nucleons carry the entire plus-momentum, 
P + of the nucleus |)15|| . This result is obtained independently of the specific relativistic mean 
field theory used, so no such theory contains the binding effect discussed above. The only 
binding effect arises from the average binding energy of the nucleus (16 MeV for infinite 
nuclear matter), and is far too small to explain the observed depletion of the structure 
function. This is in accord with an earlier similar finding by Birse [JIB]] . The generality of 



this result encourages us to seek a broader context. This is found in the Hugenholtz-van 



Hove theorem |17j which states that the binding energy of the level at the Fermi surface is 
equal to the average binding energy, or the energy of the level at the Fermi surface Ep is 
equal to the nuclear mass divided by A: 

E F = M A /A = M. (1.2) 

This theorem is the consequence of using the condition that the total pressure of the nucleus 
vanishes at equilibrium, and the assumption that nucleons are the only degrees of freedom 
contributing to the nuclear energy. Thus this theorem is a signal that P + = P^ or that 
nucleons account for the entire plus momentum of the nucleus. This generally is understood 
to imply that there will be no EMC binding effect || , thus any model which obeys Eq. ( |1.2| ) 
can be expected not to have one. 

The next step is to recall in Section |TTT| how light front dynamics is applied to computing 
the properties of infinite nuclear matter using the Walecka model (as a specific example) in 
mean field approximation (MFA). The purpose is to illustrate the general formalism needed 
to go beyond the mean field approximation, provide an explicit example of the general 
results presented in Section II, study the nuclear structure origins of the nuclear momentum 
content, and show explicitly that the Hugenholtz-van Hove theorem is satisfied. 

In Section [IV| we introduce four other model Lagrangians in which the values of the 

effective mass and vector meson field vary widely. Again our specific calculations are limited 

to the MFA. However, in Section |V|, the application of the Hugenholtz-van Hove theorem 

i7| allows us to make some general statements about models which include nucleon-nucleon 



correlations. In particular, we use this theorem to explain why no binding effect is contained 
in any model, such as that of Ref. |]18||, in which nucleons are responsible for the entire plus 



momentum of the nucleus. This is in accord with early observations of Ref. ||, but now 
there is an additional ability to compute all of the relevant nuclear properties using light 
front dynamics. We also use our findings assess existing treatments of the binding effect. 
Section [VI| is a summary of our results and their implications. In Appendix |A] we use light 
front dynamics to compute nuclear properties of the four models of Section [IV]. 




FIG. 1. Diagram for computing the nuclear structure function. A nucleus of momentum P 
emits a nucleon of momentum k, which emits a quark of momentum p, which absorbs the virtual 
photon of momentum q. 

II. DEEP INELASTIC SCATTERING FROM NUCLEI 



We are testing the hypothesis that conventional nuclear dynamics can explain the EMC 
effect. This means that we need to include possible binding energy and Fermi motion effects, 
but not dynamics related to true modifications of the nucleon structure or off-shell effects 
caused by the nuclear medium. The key assumption is that the system formed by the 
absorption of the photon is not a bound nucleon and therefore does not have the same 
interaction. The relevant lifetime of the struck system is —7? < 0.5 fm (for x > 0.5) which 
corresponds to a very short nuclear time, too short for interactions. In this case the use of a 
manifestly covariant formulation to derive the expression for the structure function leads to 
a convolution formula. If one uses the free nucleon structure function (neglecting off-shell 
effects) and F 2 n for free nucleons one finds |19 | 
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where P is the total four momentum of the nucleus, and 



x A = Q 2 A/2P ■ q = xAMjM A = xM/M 



(2.1) 
(2.2) 

(2.3) 



with M as the free nucleon mass and M = Ma/ A. The function x(fc, P) is (proportional 
to) the connected part of the nuclear expectation value of the nucleon Green's function 
p0|| , and the trace is over the Dirac and isospin indices. That /jv(y) is a Lorentz scalar 
is manifest from the structure of Eq. (|2.2|). We note the appearance of M instead of the 
free value of the nucleon mass M. This arises here from the definition ( |2.3|) and the feature 
that, in the Bjorken limit, the nuclear structure function depends on the ratio p + / P + = 
(p + /k + )(k + / P + )i where p^ is the quark momentum, with P + = Ma = AM. The basis of 
the formula ( |2.2j ) is that both the quark and nucleon distributions are directly related to 
manifestly covariant Green's functions [|19| , [2Tf . This is a standard result using nothing more 
than the stated assumptions and the Feynman diagram in Fig. 1. 



The manifestly covariant form of the single nucleon Green's function has been known for 
a long time p2 |, and its use (in the nucleus rest frame) leads to the result 
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The general form of the Green's function depends on a vector potential V = {V°, 0) for a 
nucleus at rest, and the effective mass M* which includes the effects of interactions on the 
nucleon mass. The values of V and M* depend on the specific Lagrangian employed, but the 
form of the Green's function is general. Recall also that V~ = V + = V° for the expectation 
values of vector meson fields in the nucleus rest frame. 

The result ( |2.4| ) was first obtained using the conventional equal time approach, but the 
very same can also be obtained from the light front formalism. In that case it is necessary 
to include the effects of the instantaneous part of the nucleon Green's function and those of 
the instantaneous meson exchange. 

The next step is to insert the connected part (second term) of ( |2.4| ) into Eq. ( |2.2| ) for 
/jv(2/)- This gives, after taking the trace and using the delta function to integrate over k°, 
the result 
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d 2 k< dk 3 — K =^. — 5{y v ; %1 )6(k F -\k\). (2.6) 
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The integration is simplified by using the transformation 

k + = E*{k) + k 3 , 



(2.7) 



which makes a connection with light front variables [2~3|. It is an exercise in geometry to 
show that the Fermi volume can be re-expressed in terms of k + using 



k\ + (k + - E F ) 2 < k 2 F , E F = sjkl + M* 2 , 
so that Eq. ( p.6|) becomes 

4 r .0. /".. ,„/.o .0 ,. , „_o\ „, k + + g v V A 
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I d 2 k ± I dk + 9 (k 2 F -k\- {k + - E* F f) 5{y - 



(2tt)3 Pb 7 " -J - ■ V-* - - -*> J'-' M 

The use of the definition of the energy of a nucleon at the Fermi surface, 

E F = E F + g v V + = E F + g v V°, 
allows one to achieve a simple expression for /jv(y): 
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(2.8) 



(2.9) 
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(2.11) 



The result Eq. (|2.11|) can be further simplified by using the Hugenholtz-van Hove theorem 
displayed in Eq. (|1.2|) . Section |(| contains an explicit demonstration of Eq. (|1.2j ) for the 
Walecka model and the Appendix contains a similar demonstration for the other relativistic 
models evaluated using the mean field approximation. Using Eq. ( |1.2|) in Eq. Q2.ll ) therefore 
leads to the general result 
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correct for any relativistic mean field theory of infinite nuclear matter. Different theories 
with the same binding energy and Fermi momentum may have very different scalar and 
vector potentials, but must have the same Jn^V)- 

A result very similar to Eq. Q2.12D was previously obtained by Birse [|j|. The difference 
between his formula and ours is the appearance of M in the function /jv(y), whereas he uses 
M. This difference is a small effect numerically, and therefore our conclusions will be the 
same as his. 

The baryon sum rule and momentum sum rules are derived by taking the first two 
moments of /jv(z/)- This gives: 



dyf N (y) = 1 (2.13) 

dyyf N (y) = l. (2.14) 

The latter equation is remarkable; it states that in deep inelastic scattering the nucleons 
act as if they carry all of the P + of the nucleus even though the mesonic fields are very 
prominent. 

This is clearer if we re- interpret these sum rules in terms of a probability /tv(A; + ) that a 
nucleon has a plus momentum k + = yM, with ,/jv(fc + ) = Af N (yM)/M, so that 

J dk + f N (k + ) = A, (2.15) 

dk + k + f N (k + ) = AM = M A (2.16) 



/ 



The momentum sum rule ( 2.16Q shows the total plus momentum carried by the nucleons (as 



seen in deep inelastic scattering) is also the total momentum carried by the nucleus. 

The main result of this is that the nuclear structure function is given by Eq. (|2.2|) with 
the function fniv) obtained in Eq. (|2.12|) . This tells us that, despite the fact that there is 
considerable binding energy, there is no EMC binding effect. Indeed, F 2 a depends on the 
Fermi momentum but does not depend on the effective mass M*. 

The quantity measured in deep inelastic scattering is the ratio defined by 

A numerical study of this expression using, five different relativistic models is presented 
below. First, we emphasize the qualitative features. Since the width of /at is given by 
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TABLE I. Summary of the Models- Taking the sum of g v V + /M and E F /M shows that each 
model satisfies the Hugenholtz-van Hove theorem, Eq. Ql.2|). 



Model g v V + g v V + /M E* F /M M*/M k F (fm^ 1 ) Pn/ p a 



w 


9vPB/ml 


0.355 


0.645 


0.56 


1.42 


0.65 


NVW 


1Gp B 


0.355 


0.645 


0.56 


1.42 


1 


ZM 


alpB/ml 


0.079 


0.921 


0.85 


1.42 


0.92 


NVZM 


1Gp B 


0.079 


0.921 


0.85 


1.42 


1 


RF 


ZGrfpb 


0.194 


0.806 


0.73 


1.31 


1 



the small quantity kpjM it is a very narrow function. In this case, one may evaluate the 
integrand of Eq. ( |2.2|) by expanding F 2N (x/y) in a Taylor series about y = 1 |23j to find 



that 

R{x) = Fh^ + kl (2 XaF > n ( Xa ) + x A F>' N (x A )) , (2.18) 

F2n{x) ioM F 2N (x) v J 

which shows that the only effect of the binding energy occurs in the small difference between 
xa and x which depends only on the small average binding energy. Note that a term pro- 
portional to F^ix) (but not proportional to the small parameter =?)) vanishes because one 



is expanding about y = 1 and using the baryon and momentum sum rules Eqs. ( p . 1 5| , p7T6| ) . 



We may further approximate R(x) by expanding the first term about the value xa = x 
(xa = 1-02 x for nuclear matter, and xa < l.Olx for finite nuclei). Thus 

R W = l + WjFr\ + ,,-Jl ; ; (pAFUixA) + xIF^xa)) , (2.19) 

M F 2N (x) io M F 2N (x) v ' 

where (e) is the binding energy per nucleon (16 MeV for infinite nuclear matter and < 8 
MeV for finite nuclei). This shows that, as long as the Hugenholtz-van Hove theorem is 
applicable, the only binding effect is due to the binding energy per nucleon. One is not 
allowed to use the separation energy which is much larger. The use of Eq. ( p,19|) cannot 



lead to a large enough depletion of R(x) [|l(| to resemble the extrapolated data for nuclear 



matter 25 



The relevant parameters of these models are displayed in Table |I[ The qualitative 
features discussed above are prominent in the numerical calculations displayed in Fig. |2| 
in which the ratio R(x) of Eq. (|2.17|) is presented for five different relativistic models. The 
relevant parameters of these models are displayed in Table |. We note that four of the models 
have identical results, with only the Rusnak-Furnstahl model (with its chosen different 
value of kp) differing only very slightly. These results are obtained using a simple early 
parameterization [Q of the free nucleon structure function 

F 2N (x) = 0.58^(1 - x) 2 - 8 + 0.33^(1 - x) 38 + 0.49(1 - x) 8 , (2.20) 

but the essential features of the curves are independent of the free nucleon structure function. 
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FIG. 2. R(x) vs. x for W, NVW, ZM, NVZM (solid line) and RF (dashed line). The data 
shown here are from the extrapolation of Ref. [25] 



III. THE CANONICAL FORMALISM WALECKA MODEL 

We illustrate the canonical light front formalism using the Walecka model as our first 
example. Some of this has been published previously ||, but our purposes here are to set up 
and illustrate the formalism necessary to go beyond the mean field approximation, provide 
an explicit example of the general results presented in Section II, study the nuclear structure 
origins of the nuclear momentum content, and show explicitly that the Hugenholtz-van Hove 
theorem is satisfied. 

The Walecka model employs a Lagrangian containing fields for nucleons (ip'), scalar 
mesons <fi and vector mesons V 1 : 

C w = l -{d^d,<P - m 2 s <p 2 ) - ~^%„ + \mlV^ + ^ (y(ify - g v V,) -M- g s cj>) rj;', (3.1) 



with the field equations: 



(<V + m 2 s )<P = -grfif/ 



The symmetric canonical energy-momentum tensor is given by |]2T|-f29 



(3.2) 
(3.3) 
(3.4) 



T v» = _^»£ + V cvy(}v gj3a + Qn^ffSQ + - 1 //('f(iff / - g v V u ) + Y{^ - g v V»))ij'. (3.5) 



A. Mean Field Approximation for Infinite Nuclear Matter 



We follow the MFA |T0[ in assuming that the sources are sufficiently strong so that 
the resulting large numbers of mesons can be treated in a classical manner in which source 



operators are replaced by their expectation values in the nuclear ground state. Furthermore, 
for a system of infinite volume, all positions and directions (in the nuclear rest frame) are 
equivalent. In this case <fi and V° are constants and V = 0. The approximate mesonic 
equations of motion become 

<p = -— 2 m') = -— 2 Ps (3.6) 

v° = v± = 4<v'V> = ^pb, (3.7) 

in which the brackets are used as an abbreviation for taking the ground state matrix element, 
and 

p B = 2k 3 F /3ir 2 . (3.8) 

The nucleon, though described in terms of four-component spinors, consists of only two 
independent fields. The independent and dependent degrees of freedom are defined by the 
projection operators: A± = §7°7 ± , ip' ± = A±ip', with ip' + chosen as the independent field. 
One more step is necessary because the resulting equation for ip'_ depends on V + in a 
complicated manner. It is traditional in light front dynamics of mass-less vector bosons to 
remove the effects of the term V + , by working in a gauge with V + = 0. Here we use the 
Soper-Yan transformation [27LB5|: 



tf = e~ i9 ^, d + A = V + , (3.9) 

which allows a simple equation for if)'_ in terms of ip' + to proceed in a satisfactory manner, 
but which also causes the loss of manifest covariance. With this transformation the final 
version of the nucleon field equation becomes 

{id~ - g v V-)^ + = (a ± ■ p ± + [3(M + g s (j)))^- 

id + ^ = (a ± -p ± + (3(M + g s( j)))iP + (3.10) 

within the mean field approximation (in which d~A = 0) for infinite nuclear matter. The 
nucleon mode functions are plane waves so ip ~ e tk ' x and 

{id- - g v V-)^ + = X + k+ 1>+, (3.11) 

where M* = M + g 8 <f>. 

The relevant components of the energy-momentum tensor in the mean field approxima- 
tion MFA are obtained by using the constant meson fields of Eq. ( |3.6j ) and ( |3.7| ) in Eq. ( |3.5| ) 
to obtain 

Tti + FA = m 2 X + 2^«9 + ^+ (3.12) 

TtfFA = ™ 2 S <P 2 + 2i/>l{id- - g v V-)^ + (3.13) 

and so P + and P~ are given by 
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P ± = (T& A )Sl, (3.14) 

where Q is the volume of the system, taken as infinite at the end of the calculation (A, Q — > oo 
with A/Q finite). The evaluation of these expectation values yields 

r = ^ + ^/,^* f ^F Q! (3 ' 15) 

^ + =mX + j^- 3 J F d 2 k ± dk + k + . (3.16) 

It is necessary to define the Fermi sea within the present context. Although we do not have 
manifest rotational invariance here, this invariance is restored in the results if we define the 
component k 3 implicitly through Eq. (|2.7| ). Then 



J d 3 k--- = f d 3 k6{k F - k) ■ ■ ■ , (3.17) 

and geometry leads to 

f d 2 k ± dk + ■■■ = f d 2 k ± dk + 6 (4 -k\- (k + - E* F ) 2 ) ■■■. (3.18) 



Using Eqs. ( |3.15| - |3~T7|) leads to the results that the value of the energy of the system in the 
rest-frame, 

E A = ±(P + + P-), (3.19) 

is the same as in the usual treatment of the Walecka model, as shown below. The only 
remaining task is to determine the Fermi momentum, kp. This is done by using the mini- 
mization 

E A (k F ) = M A . (3.21) 

Carrying out the differentiation leads to an equation which is equivalent to setting P + = 
P~, which must occur for a system in its rest frame with P 3 = 0. Since rotational invariance 
is maintained in the solution P 1,2 > 3 = 0, and therefore the pressure P = 1/3 J2i=i,3 P % 
vanishes. Thus the equation P + = P~ is also the light front equivalent of setting the 
pressure P to 0. Note also that one may explicitly carry out the differentiation to find that 

E A /A = M A /A = E* F + g v V° = E F (3.22) 



which is the Hugenholtz-van Hove theorem fll~T 



The above paragraph is serves as an outline of the derivation of the Hugenholtz-van Hove 
theorem, but we also provide an explicit proof. First, use the transformation ( |2.7| ) to obtain 
the results 
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s~= m «* 2+ (2^i> (£ *«4 fc - fe > (3 - 23) 

i4^ 2+ ^^ + (4// tew <3 - 25) 

Next carry out the differentiation in Eq. ( |3.20| ), using A = ps$l to obtain 

okf k F 

The term Qjjl£ is obtained by first eliminating all derivatives with respect to <\> using the 
feature that setting -J^ 5 - to zero reproduces the field equation for <fi. Then one uses the field 
equation for the vector meson (|3.7f ). The result is 

_1_|,3 B . = ^ _ ^ + _1_ / p A B . W , (3.27, 

This is a transcendental equation which determines kp, so that the calculation of E& is 
complete. With the self-consistent Eqs. (|3.6| , ^7| , ^.25|) one obtains an average binding energy 



of 15.75 MeV with a Fermi momentum of kp = 1.42 fm 1 using the parameters: -%M 2 = 
195.9, ^M 2 = 267.1, which corresponds to g v V~ = 323 MeV, and M*/M = 0.56. These 

Trig 

parameters are the same as in the original Walecka model. 

The relation P + = P~ (which must hold for a system in its rest frame) also emerges as 
a result of this minimization. To see this, rewrite the left hand side of Eq. (|3.27| ) as 



4 Air 



^-(2^// 3 *H + 3^))- < 3 - 28 > 



(2vr) 3 3 
Using this in Eq. ( |3.27| ) leads to 

m 2 , (0 rn 2 r0 1 r ,, , k ■ k 



l^v 2 = ——r / d 3 k '" ,'" , (3.29) 

2 ° (2tt) 3 Jf 3E*(kY V ; 

which is what one obtains by setting P + = P~ using Eqs. ( |3.23| ) and ( |3.24j ). 

We now have the tools at hand to prove the Hugenholtz-van Hove theorem. Simply use 
Eq. ( |3.27D to remove the integral appearing in Eq. ( |3.25| ) and obtain 

^ = m 2 v V 2 + p B E F . (3.30) 

Then the use of the field equation (|3.7| ) yields 

h ' ^=g v V° + E F = E F , (3.31) 



p B fi A 

which is the desired result. This is a remarkable result. The original version of the theorem 
was proved using only the assumption that nucleons are the only degrees of freedom. Here, 



the mesons are important, yet the theorem still holds [[!(J 
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B. Nuclear Plus-Momentum Content 

Now we relate the role of the plus component of the momentum seen here with that of 
Section |I[The nucleonic contribution to the nuclear plus momentum from Eq. ( |3.16| ) is 

p+ 
N - * d 2 k ± dk + k + , (3.32) 



A pb(27t) 3 Jf 

which is also obtainable directly from taking the nuclear expectation value of nucleon plus 
momentum operator. 

The large vector potential and small value of M* are associated with the startling result 
that only 65% of the plus momentum of the nucleus is carried by nucleons, and that 35% 
is carried by vector mesons. It was previously argued |7||| that this would produce a dis- 
astrously large decrease in the nuclear deep inelastic structure function. As shown above, 
that does not occur, because the function fn{y) peaks at y — 1. 

We therefore need to understand how it is that the nucleons can carry 65% of the 
momentum here all of the momentum as stated in Section [TI|. To do this, use Eq. ( |3.32j ) to 
define a probability f{k + ) that a nucleon carries a plus momentum k + : 

P+/A = f dk + k + f{k + ), (3.33) 



with 



f ^ = —FTVs I ^ = —FTVs I d2k ^ d P + ^ k+ ~ P + )- ( 3 - 34 ) 

Pb{2ix) 6 Jf pb(2it) Jf 



Pb(2tt) 3 Jf pb(2tt) 

It is useful to again obtain a dimensionless distribution function f(y) by replacing k + by 

fe+ 

M 



the dimensionless variable y using y = =, f(y) = Mf{k + ). Then one finds 



3M 3 



f(y) = l—0(y + -y)9(y- y - 



4 k F 



k 2 R* 

M 2 { M V) 



(3.35) 



where y^ = F — F . This function peaks at y — E* F / M = 0.65 for the Walecka model, and the 
average value of y is also 0.65. Its use in Eq. ([2.1|) would indeed lead to a disaster. Indeed, 



our previous work assumed that Eq. (|1.1|) was appropriate. In that case, the computed 
ratio 2A ^ X ' was dramatically smaller than F 2 n- This disastrous result can be understood 
from the following logic. A reasonable first approximation to the integral Eq. ( |1 . 1| ) can be 
obtained by using f(y) w 5(y — ■=■) which satisfies the baryon sum rule and corresponds 

to an average value of y = 0.65. Then 2A ^ X ' vanishes for x > 0.65 and the ratio to the 
free structure function goes to zero in huge contradiction with experiment, which shows 
depletions no larger than 20% for the heaviest nuclear targets. This result is illustrated with 
the numerical calculation shown in Fig. |3|. 

However, the correct quantity to use in deep inelastic scattering is /v(y), which emerges 
from a manifestly covariant treatment. To see the connection between f(y) and /jv(y), it is 
first helpful to compare directly Eqs. ( |2.12|) and ( |3.35| ). This comparison yields: 
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FIG. 3. R(x) vs. x for the Walecka model using fN(y) (solid line) and f(y) in Eq. ( |2.lD (dashed 
line). The data shown here are from the extrapolation of Ref. pH] 



f(y) = fN(y + g v v + /M), (3.36) 

using the Hugenholtz-van Hove theorem Eq. ( |1.2| ). The correct nucleon distribution function 
is one that is shifted by the vector potential. This simple relation suggests that there is a 
simple interpretation of the difference between JN{y) and f(y) in terms of a phase difference. 
Indeed, the difference arises from the Soper-Yan transformation ( |3.9|) which relates the fields 



tj) and ip' . Thus we have two forms of the plus momentum density, T 1 



++ . 

MFA- 



1 MFA 



m 2 X + 2^ + id + ij + = 2i/>'* + (id + - g v V + ) ^ + 



(3.37) 



In the second form, a nucleon operator carries all of the plus momentum, the correct single 
nucleon plus momentum operator is the canonical conjugate momentum which is shifted 
by the term g v V + . The second form is the appropriate one as it is related to the original 
covariant Lagrangian. The term yfN{y) is obtained from the expectation value of the nucleon 
plus-momentum operator ip+(id + — g v V + )ip' + , while yf(y) is obtained from ip + id + ip + with 
/jv(y) as the distribution function which emerges from a covariant treatment. 

IV. FOUR MORE MODELS IN MEAN FIELD APPROXIMATION 

The Walecka model, evaluated in MFA, was known to have some phenomenological 
troubles. The compressibility is too large, and the very small effective mass (shown here to 
be irrelevant for deep inelastic scattering) does enter into quasi-elastic scattering, in which 
it is a straightforward matter to show that the cross section is given by an integral over the 
distribution, f(y), and not /jv(y)- The reason for this is that in the MFA the struck nucleon 
feels the same vector and scalar potentials as a bound nucleon. Hence it was of interest to 
improve the Lagrangian. This has been done in a variety of ways. We consider four other 
models here, mainly to show that the same function /j\r(y) emerges from each one and that 
each one satisfies the Hugenholtz-van Hove theorem. The validity of this theorem is a signal 
that the nucleons carry all of the plus-momentum, just as for the Walecka model, so that 
there can be no significant binding effect. 
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All of the models have essentially the same saturation properties, but each is distin- 
guished by using a different mechanism to reduce the putative amount of momentum carried 
by the vector mesons. Here we simply define the models and summarize the results. The 
details of the solution are presented in the Appendix. The important parameters of each 
model are displayed in Table |I[ 

The No Vector Walecka (NVW) model is defined by the elimination of the vector meson 
field in favor of a point coupling interaction of the form Gj^j^ with j^ = ip'j^ip'. In this 



case the nucleons carry all of the momentum. The Zimanyi-Moszkowski (ZM) model JTJ]] is 
defined by using a rescaled derivative coupling interaction in which the scalar coupling in the 
Lagrangian is given by the term — ip M/(l — a ^-)ip' '■ This model is known to involve a larger 
effective mass and smaller vector field than the Walecka model. The Appendix shows that, 
in this model, the nucleons carry about 92% of the total plus momentum. The No Vector 
Zimanyi-Moszkowski (NVZM) model is defined by starting with the Zimanyi-Moszkowski 
Lagrangian and then removing the vector mesons in favor of a current- current interaction as 
in the NVW model. Again nucleons carry all of the plus-momentum. The Rusnak-Furnstahl 
(RF) point coupling model contains no explicit meson fields, and the interactions included 



via a variety of non-linear couplings. The parameters are given in Ref . |12| , and arguments 



for their naturalness in terms of effective field theory have also been presented |I3" |. The 
nucleons carry all of the plus-momentum and the numerical value of the effective mass is 
M* = 0.73M, significantly higher than that of the Walecka model. 

Each of the models summarized in Table | has essentially the same saturation properties 
even though the values of M* and E* F display huge variations. Note that each model has 



a nucleon mode equation ( |3.11| ), (|A9|) , ( |A16| ), ( |A29| ) , and ( |A36| ) , and these are summarized 



by the single unifying expression: 

(id- -g v V+)ip + = kl+ ^* ] \ + . (4.1) 

The numerical values of g v V + are listed along with other relevant information regarding the 
five models in Table |. 

V. BEYOND THE MEAN FIELD APPROXIMATION 

It is worthwhile to discuss the generality of the result that there is no significant binding 



effect. Consider any model, such that of Ref. [|18j in which mesonic fields are not explicit 
components of the nuclear Fock state wave function. For example, one may eliminate the 
mesonic degrees of freedom in favor of two- and three- nucleon interactions without main- 
taining the mesonic presence in the nuclear Fock state wave function. Such models, correctly 
evaluated, obey the Hugenholtz-van Hove theorem, Eq. ( |1.2| ). The validity of this theorem 
is a signal that nucleons carry all of the plus momentum, so that the baryon and momen- 
tum sum rules Eqs. ( |2.13| - |2~7L6| ) are satisfied. This means that for equilibrium the following 
conditions hold 

P ± = P^ (5-1) 

P+ = P~. (5.2) 
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These two equations (the first is defined by the model, the second by stability) may be 
thought of as the light front version of the Hugenholtz-van Hove theorem. Therefore one 
may again apply the analysis of Refs. [0,0 and expand /jv(y), appearing in the integral of 
Eq. (|2.2| ) about its peak value of unity It is generally sufficient to keep only three terms. 
Thus one finds 

F 2A {x) = F 2N {x A ) + 7 (2x A F^ N (x A ) + x 2 A F^ N {x A )) (5.3) 

7= fdy(y-l) 2 f N (y). (5.4) 



The coefficient 7 is larger than the term proportion to k 2 F of Eq. ( |2.19| ) because the effects 
of correlations extend the width of the distribution /jv(?/), but it multiplies a term which is 
positive in the valence quark region. Thus, once again we see that the only binding effect 
appears in the presence of the variable x A which is only slightly larger than x, see Eq. ( [2.3| . 
This effect is too small to reproduce the data; there is essentially no EMC binding effect. 
The result (|5.3|) is very similar to Eq. (|2.19|) in that the term, usually associated with the 
binding effect, proportional to ■§-F 2 N{x A /y)\ y =i vanishes because of the second sum rule of 
Eq. (|1.2j ). We stress that any model in which nucleons are the only degrees of freedom must 
obey the momentum sum rule, as expressed in either Eq. ( [2.14| ) or fl2.16| ). Thus Eq. ( p.3| ) 



will emerge and the model will not have a sufficiently large binding effect to explain the 
nuclear deep inelastic scattering data at large values of x. 

Can the conventional meson-nucleon picture of nuclear structure (which ignores off-shell 
effects) be used to reproduce the nuclear deep inelastic scattering data in the valence quark 
region of Bjorken xl The only way to get a binding effect is to compute the nuclear ground 
state wave function in such a manner as to obtain the mesonic P^ and nucleonic contribu- 
tions: 

P + = P+ + P+, (5.5) 

P- = P N +P m , (5.6) 

in which the meson content P^ is treated explicitly. In general, P^ consists of terms arising 
from any of the exchanged mesons which are responsible for the nuclear force: 

P* = P* + P± + P* : + ---. (5-7) 

The equation for P + was used long ago [0^1 , in which nucleons and pions contributed to 



the total, and with momentum conservation presented as the justification for the equation. 
It is useful to realize that the use of the energy-momentum tensor provides a general basis 



for this sum rule. For example, in the work of Refs. ||3|,|3~Tf the use of a chiral Lagrangian, 
containing isoscalar vector mesons, V* 1 , V^ u = d^V u — d u V' 1 , scalar mesons <fi and pions 7r, 
and standard manipulations give 

T++ = v ik V ik + m 2 v V + V + + ^ + id + i) 

+ d + <f)d + <P + d + 7r -d+ir + ir- (9+tt 2 ^ (1 - ^sin 2 ^), (5.8) 



IT 2 7T 2 /' 



with P + given by 
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P+ = (T ++ )Q, (5.9) 

with the brackets denoting a ground state matrix element. The point is that each separate 
term corresponds to a term in either P^ or P+, and therefore identifiable as a term in 
the sum rule fl5.5|) . However, the field equations provide relations between all of the fields. 
Thus, one can not get a reasonable result for P + by considering only one of the terms which 
contribute. 

The pressure balance condition P + = P~ = AM, must hold for a stable solution so that 
one finds 

P N — P N = P m — P m . (5.10) 

Thus the condition needed to prove the the Hugenholtz-van Hove theorem (that P^ = P^ = 
AM) is not obtained. We know that P£ < Ma because all of the contributions to P+ are 
positive definite. One may therefore define a positive quantity e via the deviation: 

es S- (5 ' n) 

so that 

Jdyyf N (y) = l-e. (5.12) 

Thus Eqs. ( |5.10| , [5TT|) can be thought of a generalization of the Hugenholtz-van Hove theorem 
which is equivalent to the momentum sum rule. With this new feature, the application of 



the expansion procedure to the integral of Eq. (2.1) yields a term proportional to e: 



F 2A (x) = F 2 n(x a ) + ex A F^ N (x A ) + 7 {2x A F! 2N (x A ) + x A F% N (x A )) . (5.13) 

Equation ( |5.13| ) corresponds to the usual binding effect which now is present. However, one 
needs fairly large values e ~ 0.05 to reproduce the deep inelastic data for Iron ||, and e ~ 



0.07 to reproduce the nuclear matter extrapolation shown in Fig. |2|. Early calculations ||33|J32 
in which pions are allowed to carry such a momentum fraction have another consequence [34] : 
an enhanced nuclear anti-quark content which turned out to be in contrast with the results 
of the nuclear Drell-Yan experiment [|35] . A more recent light-front calculation |36|J3^| , which 



included the effects of nucleon-nucleon correlations as well as those of an explicit meson Fock 
space, finds that pions carry about 2% of the nuclear plus momentum. While this value is 



consistent with the Drell-Yan experiment [38], using e = 0.02 would provide too small a 



reduction in the nuclear structure function. It is possible that other mesons could supply 
significant contributions to e, and it is necessary to investigate this possibility. However, 
consistency with the Drell-Yan experiment must be maintained. While it seems unlikely 
that a careful calculation will be consistent with both the deep inelastic and Drell-Yan data, 
we cannot rule that out now. 



The analysis of the previous paragraph is similar to the early one of Refs. |3j and [32 



The main difference occurs in the present ability |J to compute the nuclear binding energy 



m terms of P^ and P^. 
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It is also worthwhile commenting on the modern calculations |J which use the nucleon 
spectral function S(p) to compute the quantity fn{v)- It is necessary to obtain Eq. (|5.12|) 



with a significantly large value of e to achieve agreement with data. However, a complete 
calculation should also obtain the very same value of e from Eq. (|5.11| ). But models in 
which the mesons are eliminated in favor of two and three nucleon potentials, forfeit the 
ability to compute the value of e directly. It is possible to make a completely accurate 
calculation of fn{y) which would reproduce the correct value of e, but a computational 
error which is only a few percent in /jv(y) corresponds to a huge percentage error in the 
small quantity e. Hence, such calculations must be regarded as inconclusive. Even if we 
take the calculations at face value, the models "are not completely satisfactory" . If mean 
field models are used, nuclear binding accounts for only 20% of the observed effect |39||.Very 
large separation energies (values of e), inconsistent with the mean field calculations, are 



required PUB to reproduce the data. Calculations have been made including correlations, 
but the summary of Refs. [0,0 made in the book || is: "But for all nuclei considered 
the predicted deviation of the ratio R(x) (is) much smaller than the experimental one." 
This statement applies also to the work of P5| , if the "off-shell nature of the nucleon" 
is ignored. The inclusion of off-shell effects by allowing the nucleon structure function to 
depend on the momentum of the nucleon in the nucleus (as well as on x/y) can lead to 
a significantly improved description of the data [H3,H3|. This agreement is consistent with 



the results of the present work. Here we consistently ignore off-shell effects, categorizing 
these, along with a host of others, as interesting effects. In any case, one would need to 
understand the implications of analogous off-shell variations in operators used in impulse 
approximation calculations for many nuclear reactions. Furthermore, it is not clear to us 
that these formulations |4^j43|] provide nuclear structure functions which are consistent with 
the baryon sum rule. 



VI. SUMMARY AND DISCUSSION 

The principal result is that relativistic mean field models of nuclei, successful for many 
observables, do not contain the binding effect needed to reproduce the depletion observed 
by the EMC. The generality of this conclusion is related to the use of the mean field approx- 



imation, consistent with the Hugenholtz-van Hove theorem |T7|], which severely constrains 
the nucleon distribution function /jv(y). This theorem has the further implication that any 
model in which the entire plus momentum is carried by the nucleons, in the sense of Eq. ( |1.2[ ), 
also contains no binding effect. Thus including nucleon- nucleon correlations (within a model 
containing only nucleons in the Hamiltonian) cannot reproduce of the data. A minimum 
feature necessary to describe the data using conventional meson-nucleon dynamics is that 
the mesonic components must comprise an explicit part of the nuclear Fock state wave func- 
tion, and the mesons must carry a significant fraction of the nuclear P + . But there are 



severe constraints on the nuclear anti-quark content @,@ and these limit the flexibility of 
mesonic models. Therefore, all of our present considerations are consistent with the notion 
that some effect not contained within the conventional framework is responsible for the EMC 
effect. 
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APPENDIX A: FOUR MORE MODELS 

Four other relativistic models are solved using light front dynamics in this section. The 
techniques are the same as in Section [HI] , so our treatment will be briefer than what appears 
above. 

1. The No Vector Walecka Model 

To reduce the effects of vector mesons it is natural to employ a Lagrangian that has 
only one meson field, a scalar field 0. The repulsion is supplied by a repulsive vector point 
coupling. The Lagrangian is 

C WNV = \{d^d^ - m 2 ^ 2 ) + W{ip~ M - g s <j>W - Gf j„ (Al) 

where 

U = W-yrf. (A2) 

The resulting equation of motion for the Dirac field is 

y(id„ - 2Gj^' = (M + gJW, (A3) 



and the equation for the scalar field is again Eq. ( |3.2| ) The canonical energy-momentum 
tensor is given by 

T*" = -g^C + d^(j)d v (j) + ^jP{-fer + -fVW. (A4) 

In the Mean Field Approximation (MFA) the equations of motion are given by Eq. ( p.6| ) 
for the scalar field and 

^(id, - 2G{ 3 ,)W = (M + g 8 W, (A5) 

for the Dirac field. The components of the energy-momentum tensor are given by 

TtiU = iWl + d + i>' (A6) 

TtfFA = m 2 J 2 - 2 Wirfiid* ~ G(jJ) -M- g^W + ^W(7 + ^ + i~d + W (A7) 

We solve the Dirac equation as in Section [Tfl| using the transformation 



V>' = e" 2iGA(a; V, <9+A = (j + > (A8) 

and define j M = (j M ) — <9^A (note that j + = by construction and j J = in the rest frame 
so that the only non- vanishing component in the MFA is j~ = ps). The Dirac equations for 
ip+ and i\)_ become 

(id~ - 2G]-)^ + = (a ± ■ p ± + {3(M + g^))^ 
id+^ = (a± • p± + (3(M + g s (j)))ip + 

If we assume ip ~ e lk ' x , then we obtain 

en- on~-\i kl + {M + g s 4>) 2 

(id - 2Gj )ij) + = — 1/>+. (A9) 

Returning to the energy momentum tensor, which has also changed under the transformation 
Eq. (|58D , we find 

T++ FA = 2^ + (id + + 2G(j + »^+ (A10) 

TtfFA = m 2 J 2 + 2^ + (id~ - 2G]-)^ + . (All) 

Using Eq. (|A9|) in Eq. ( |A11| ) and the light front 4-momentum definition Eq. Q3.14J ) we 
obtain 

p- 

,2/2 



Q i2n) 3j F d2k ^ dk+ ( k+ + 2G PB) (A13) 



— = m 
Q 

P + 



The second term in Eq. ( |A13| ) can be rewritten 

4 f J2, „+o^ 8G Pi! 



/ d 2 k ± dk + 2Gp B = —4? / d 2 k ± dk + = 2Gp] 
Jf (27rr Jf 



(27t) 3 Jf " rD (2 

and the resulting equations are exactly those of the Walecka model and we draw the corre- 
spondence (with kp = 1.42 fm _1 ) 

n 2 
2GM 2 -► ^-M 2 , (A14) 

m 2 

which means that the saturation properties of this model are the same as those of the 
Walecka model. 

An important difference between this model and the Walecka model is that the extra 
term is not due to vector mesons, but part of the nucleon contribution to P + . All of the 
plus momentum is due to the nucleons and not the vector mesons. It is apparent that this 
model is consistent with the Hugenholtz-van Hove theorem. The values of P ± are the same 
as those of the Walecka model for all values of k F . 
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2. The Zimanyi-Moszkowski Model 

The rescaled derivative coupling Lagrangian given by Zimanyi and Moszkowski J51J is 
known to have a smaller vector potential than the Walecka model. This model is defined by 
the Lagrangian 

Czm = \{d^d,<P - m 2 2 ) - \v^V, u + \mlV»V, 

+^ (V(z<9 M - g v V,) - -^\ 1>. (A15) 



Using the methods described in Section |T|, one finds that the eigenvalue equation corre- 
sponding to Eq. ( |3.11| ) is given by 

k 2 a. I M*\ 2 
(id' - g v V-)^ + = ± + fc l + ] ^ (A16) 

where 

M 

M* = -r (A17) 

L M 

V- = V = 9 -^-. (A18) 

ml 

The field V~ is transformed according to Eq. ( |3.9|) . The relevant components of the energy- 
momentum tensor are 

TtiU = m 2 X + 2^\td + ^ + (A19) 

TtfFA = ™ 2 s<P 2 + 2i/>Ui&- - g v V-)4> + (A20) 

and so P + and P~ are given by 

r= m ^ + (2^// 2 ^* + * + <A22) 

which are superficially the same as the Walecka model, but differ in that now M* is given 
by Eq. (|A17|) . The parameters of the model are obtained by minimizing the total energy at 
kp = 1.42 fm _1 , using the values 

a 2 

^-M 2 = 43.2 (A23) 

m l v 

a 2 

^M 2 = 140.4, (A24) 

mj 

which corresponds to 
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M* = 0.85M. (A25) 

The plus momentum may be decomposed using P + = P+ + P£ (meson part and a 
cleon part) 



p+ 

m 

A 
P + 


m 2 V 2 
__ m v v o =73MeV 

PB 
- 1 4 f rf 2 h,r1h + h + - 


= 850 MeV . 


A 


p B (2tt) 3 Jf " ' J 



M* = 0.85M. (A25) 

momentum may be < 
nucleon part) 

P+ m 2 V? 

(A26) 

(A27) 

The nucleons carry about 92% of the total plus momentum in this model. Despite this, the 
Hugenholtz-van Hove theorem is satisfied because the expressions for P^ in terms of V° and 
M * are the same as those of the Walecka model. The only difference is the relation between 
M* and <fi. However, that relation does not enter in the derivation of the Hugenholtz-van 
Hove theorem presented in Section III. 

3. A No Vector Zimanyi-Moszkowski Model 

The next step is to modify the Zimanyi-Moszkowski Lagrangian by removing the vector 



mesons in favor of a current-current interaction as in Section A 1 This Lagrangian is 



-G^H>~i^. (A28) 



The results follow exactly as those of Section [Al| . Specifically, changing from a "vector" 
to "no vector" model does not affect the minimization of the energy density and therefore 
leaves the coupling constants ( |A23| ) and ( |A24j ) unchanged if we make the identification 



2GM 2 — > (gl/m 2 )M 2 . The operator corresponding to Eq. (|A9|) is 



k 2 4- ( M*) 2 
(id- - 2G]-)^ + = x + fc l + ] ^ + (A29) 



For this model P is exactly Eq. ( |A21| ) and we have plus momentum 

P+ 4 



a, " (2tt) : 



■ f d 2 k ± dk + (k + + 2Gp B ) (A30) 



This model obeys the Hugenholtz-van Hove theorem because the values of P ± are the 
same as those of the Zimanyi-Moszkowski model for all values of &p. 

4. Rusnak-Furnstahl Point Coupling Model 



The point coupling Lagrangian of Ref. |T2| is the modern version of the original Walecka 



which is connected to QCD through symmetry and naturalness. This model is defined in 
terms of the following densities: j^ = ip'^^ip', p s = ip'ip' and s^ = ip'a^ u ip', so that 
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Cut = ip'(ip~ M)iJ/ - P 2 s (k 2 + k 3Ps + K 4 p 2 ) - j^{Q 2 + VlPs + mP 2 s + ( 2jfl f) 
-d^p s d^p s (K d + aip a ) - d^j^fiCd + a 2 p s ) - fv(d^f)s^ v . 

The parameters are given in Ref. [I^|, and arguments for their naturalness in terms of 



effective field theory have also been presented |T3 



The relevant components of the symmetric energy momentum tensor in the MFA are 
given by 

T&a = #V<9V (A31) 

T&fa = -W{iP~ M)i>' + 2 P 2 s {k 2 + n 3 p s + k aP 2 s ) 

+2j^(( 2 + ViPs + V2P 2 S + Cduf) + ^W(l + d- + Td + W (A32) 

and the Dirac equation is 

[ip- M)iP' = [ Ps (2k 2 + 3k 3Ps + Ak 4P 2 s ) 

+7mJ m (2C2 + 2 VlPs + 2 V2 p 2 s + 4&j»f) 

+J^(Vi + 2V2PsW- (A33) 

Note that the various densities are constants within the MFA, we may define 

M*=M + Ps {2k 2 + 3k 3Ps + Ak aP 2 s ) + j^fa + 2 V2 p s ) (A34) 

g rf = C 2 + vips + V2pl + 2Cd^"» ( A35 ) 

and follow Section |TJ to obtain the equation 

(id~ - 2G RF j-)^ + = L+ k \ ] V>+, (A36) 

and the momenta 

p- 

— = -2k 2P 2 s - AK 3 pl - 6k 4 Ps - 2 ViPspl ~ ^]2p\p\ - a CaP b 

^ =2G RF p 2 B + -^J F d 2 k ± dk + k + (A38) 

The derivation of Eq. (|A37| ) , involves adding and subtracting M — M* in order to use the 
Dirac operator ( |A36| ) in obtaining the last term. The numerical value of the effective mass 
is computed to be M* = 0.73M for k F = 1.31 fin -1 . 

The proof of the Hugenholtz-van Hove theorem is obtained by going through the steps 
analogous to those of Eqs. (|3.19|) - (|3.31|) for this updated model. One can see numerically 



in Table | that the model obeys the theorem. Additionally, the total energy, given by 
substituting Eqs. (|A37|) and (|A3S| ) into Eq. (|3.19|) , is due entirely to nucleons, and therefore 
the model must obey the Hugenholtz-van Hove theorem. 
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